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ABSTRACT 

We have selected pre-main-sequence stars in 12 groups of notional ages ranging from 
1 Myr to 35 Myrs, using heterogeneous membership criteria. Using these members we 
have constructed empirical isochrones in V, V-I colour magnitude diagrams (CMDs). 
This allows us to identify clearly the gap between the radiative main sequence and the 
convective pre-main-sequence (the R-C gap). We follow the evolution of this gap with 
age and show that it can be a useful age indicator for groups less than ~ 15 Myrs old. 
We also observe a reduction in absolute spreads about the sequences with age. Finally 
the empirical isochrones allow us to place the groups in order of age, independently 
of theory. The youngest groups can be collated into three sets of similar ages. The 
youngest set is the ONC, NGC6530 and IC5146 (nominally 1 Myrs); next Cep OB3b, 
NGC2362, A Ori and NGC2264 (nominally 3 Myrs); and finally a Ori and IC348 
(nominally 4-5 Myrs). This suggests Cep OB3b is younger than previously thought, 
and IC348 older. For IC348 the stellar rotation rate distribution and fraction of stars 
with discs imply a younger age than we derive. We suggest this is because of the 
absence of O-stars in this cluster, whose winds and/or ionising radiation may be an 
important factor in the removal of discs in other clusters. 

Key words: stars:evolution - stars:formation - stars: pre-main-sequence - tech- 
niques: photometric - catalogues - (stars) Hertzsprung-Russell H-R diagram 



1 INTRODUCTION 

Most ages for stellar groups, sub-groups, clusters or asso- 
ciations are derived from fitting observations to theoretical 
isochrones. It is well known that these models often deviate 
systematically from the data. For instance if the isochrone 
fits part of a sequence it may systematicall y deviate from 
the observed sequence in other sections (e.g. iBonatto et al.l 
120041 : IPinsonneault et al.ll2004l h In addition if a sequence of 
stars is fitted to an isochrone in one colour, derived param- 
eters such as the age may be different fro m those derived 
in another colour (e.g. iNavlor et al.l [2002). In this paper 
we create empirical isochrones. These have the advantage 
of necessarily fitting the entire data sequence, allowing us 
to create a relative age ladder. We can then derive ages 
for fields where different regions of a sequence are available 
by selecting appropriate fiducial sequences covering the re- 
quired colour range. The obvious disadvantage of our empir- 
ical isochrones is that they are unable to provide absolute 
age information. Although this is the first attempt to cre- 
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ate such a ladder, carried out with data that uses slightly 
different photometric systems and heterogeneous selection 
criteria, the method is clearly effective. 

As we have used a range of ages to create the empiri- 
cal isochrones we are able to examine how these isochrones 
change with age. In many fields this reveals a clear gap be- 
tween the convective lower-mass pre-main-sequence (PMS) 
and the radiative stars close to or on the main sequence 
(MS). We refer to this feature as the radiative-convective 
gap (R-C gap). 

For this study we have used the optical bands V and /. 
This enables us to minimise the effect of accretion, disc pres- 
ence and chromospheric activity on our empirical isochrones. 
Accret ion will generally affe ct the flux short-ward of the B 
band (|Gullbring et al.l 19981 ) with disc emission important 
at longer wa velengths, first b ecoming significant in the H 
and K bands (|Hartmannll 1998ft . Lastly. IStauffer et all (|2003h 
have shown that magnetically generated chromospheric ac- 
tivity in young stars primarily causes perturbations and 
scatter in the B-V and U-B colours, whereas the V, V-I 
CMD is unaffected. 

Our paper is laid out as follows. The data collection and 
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reduction are detailed in Section [2] To create the empirical 
isochrones we select members in each field using a variety 
of criteria (see Appendix lAl and Section [33} and fit a curve 
through these in a V, V-I CMD (see Section [33)1. We then 
apply a distance modulus and correct for extinction (see 
Section |4]). The reader may wish to move straight to Section 
[5] Here we present the data corrected for distance and ex- 
tinction, allowing comparisons to be made with theory and 
also between sequences from clusters, groups or associations 
of different ages. Finally we combine the sequences to form 
a fiducial relative age ladder (see Section [6]), which we then 
use to find new age estimates for six fields (see Section^. In 
Section [9] we summarise our conclusions. Appendix |X] con- 
tains a discussion of the literature membership data, while 
in Appendix|B]we discuss the distances and extinctions from 
the literature that we have used. 



2 DATA COLLECTION 

The dataset used in this paper is from both new observations 
and literature sources. 

2.1 New observations 

The images for the majority of this dataset were obtained 
in BVI with the 2.5m Isaac Newton Telescope (INT), situ- 
ated on La Palma, equipped with the four EEV CCD Wide 
Field Camera (WFC). The filters used were Sloan %' and 
Harris V and B. The datasets used were collected on the 
nights of the 27th September and the 5th of October 2004, 
when the range of seeing was 1-2" and 1-1.5" on each night 
respectively. Both nights were photometric. The observa- 
tional details are given in Table [1] Standards from lLandol j 
l| 19921 ) were observed on these photometric nights, tying our 
calibration to B, V and I c Landolt standards. 

2.2 Data reduction 

The images were debiased, and then flat fielded using me- 
dian stacks of many sky flat frames. The i band images 
were then defringed using a library fringe frame. Known 
bad pixels were flagged. The data, including the standards, 
were then reduce d follow i ng th e opti mal extraction algo - 
rithm detailed in iNavlorl (|l998r i and iNavlor et ail l|2002f ). 
We used the 2 Micron All Sky Survey (2MASS) to provide 
an astrometric solution for each field accurate to ~0.1 arc- 
sec. The reduction process provided catalogues of photome- 
try with flags and uncert ainties. The flags are explained in 
iBurningham et all (|2003T l. excepting th e flag T. These are 
objects were the profile correction (see INavlor et afl 120021 ') 
fails due to a paucity of stars in the field. Therefore aper- 
ture photometry is performed using an aperture of the same 
radius as that to which the optimally extracted fluxes are 
corrected. Since aperture photometry is noisier than opti- 
mal photometry, we only carry out this proceedure for the 
brighter stars, hence the fluxes for fainter stars in the same 
image will be optimally extracted, and flagged as having 
poor profile corrections (H). 

Photometric calibration coefficients were calculated 
from the standard star observations using a similar proce- 
dure to that detailed in lPozzo et al.l (|2003T ) . For each colour 



or magnitude we allow the following free parameters: one 
extinction coeffficient per night; one colour term for each 
CCD; one zero point for each night, with offsets with respect 
to this (which do not change from night to night) for three 
of the four CCDs. The standards used had colour ranges of 
-0.30 < V-I < 2.87 and -0.30 < V-I < 2.85, with airmass 
ranges of 1.13 — 1.32 and 1.14— 1.22, on September 27th and 
October 5th 2004 respectively. 

We estimate the combined uncertainty in profile correc- 
tion and transformation to the Landolt system by adding a 
magnitude independent uncertainty. Adding 0.02 mags in 
each band gave a reduced \ 2 of approximately one in V and 
B-V. The value of reduced \ 2 was somewhat higher in V-I 
(1.4), which we suspect is due to the poor match between 
our i filter and those used by Landolt. 

In the cases where multiple fields were observed 
(NGC2264, a Ori and h and \ P er ) the catalogues we r e com - 
bined using the method described in INavlor et al.l £2002). 
Overlapping regions were used to normalise th e photometry 
to pro vide a combined catalogue. As detailed in lNavlor et al.l 
(2002) the RMS of the overlap region is a good measure of 
the accuracy of the profile correction. This suggests a mag- 
nitude independent uncertainty of approximately 0.01 mags 
which should be added to all bands when comparing objects 
well separated on the CCDs. Using the overlap region sug- 
gest an accuracy of 1% in both V and B-V and around 2% 
in V-I. We believe the higher uncertainty in V-I is caused by 
a lack of redder standards across all the CCDs. This uncer- 
tainty is not included in the catalogues. For all subsequent 
analysis the combined catalogues have been used. 

All of the new catalogues used for this paper are freely 
available from the cluster collaboration home pag^E an d the 
CDS archive (sample can be seen as Table iBlj) . For member- 
ship selection and CMD plots we exclude stars with uncer- 
tainties greater than 0.1 mags in colour or magnitude. The 
machine readable catalogues retain all the data. 

To aid navigation through this paper an index to the 
catalogues is included as Table [2] Full catalogues are pre- 
sented for each field with new photometry in addition to 
catalogues of stars which fulfill given membership criteria. 

2.3 Literature data 

Literature observations were taken in all cases (except where 
stated) in the Johnson-Cousins photometric system using I c 
(Cousins I). This raises the question of what effect the use 
of two different I filters has on our results. The difference 
is illustrated in Figure [1] The data for Sloan i, for stars in 
Cep OB3b is from this work and the data for Cous ins I c , 
for the same stars is taken from P ozzo et al. | i|2003h . The 
symbol I from hereon will be used to represent data taken 
using a Sloan i (i') filter tied to the Landolt standards, with 
7c being data taken in Cousins I c tied to Landolt standards. 
The difference ( V-I c )-( V-I) is plotted as a function of ( V-I c ) 
in Figure [T] The data can be represented by the following 
linear fit: (V-I) = 0.979(V-I e ) + 0.009. This result shows 
that two stars with an apparent V-I of 3 mags, one cool 
and unredenned i.e. a typical Landolt standard, and another 
hotter and more heavily reddened, need not have the same 

1 http: / /www. astro. ex. ac.uk/people/timn/Catalogues/description. html 
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Field 


Filter 


Exposure time (sees) ( X 1 unless stated) 


h Per 


1 


3 and 30 




V 


1, 10 and 100 




B 


2 and 20 


X * 6T 


j 


3 and 30 




V 


1, 10 and 100 




B 


2 and 20 


NGC7160 


j 


3 and 30 




V 


1, 10 and 100 




B 


2 and 20 


NGC2264 (4 fields) 


I 


2 and 20 




V 


2 and 20 




B 


4 and 40 


IC348 


I 


2 and 20 




V 


2 and 20 




B 


4 and 40 


Cep OB3b 


I 


5(x2), 30(x2) and 300(x2) 




V 


5, 30 and 350(x4) 


cr Ori (Fields 1,2 and 4) 


I 


2 and 20 




V 


2, 10, 100 and 350(x4) 


a Ori (Field 3) 


I 


2 and 20 




V 


2, 10, 100 and 350 




13 


4 



Table 1. Exposure times for each field. 



Region 


Table Number 


Data 


Source (as table [3} 


NGC2547 


4 


Members 


1 (X-ray) 




5 


Members 


2 (Spectroscopy) 


NGC7160 


6 


Full catalogue 






7 


Members 


3 (Extinction) 


NGC2264 


8 


Full catalogue 






9 


Members 


6 (Periodic Variables) 




10 


Members 


4 (X-ray) 




11 


Members 


5 (Ho) 




12 


Members 


26 (Proper Motion) 


Cep OB3b 


13 


Full catalogue 






11 


Members 


7 (X-ray) 




15 


Members 


8 (Spectroscopy) 




16 


Members 


18 (Periodic Variables) 




17 


Members 


19 (Ha) 




18 


Members 


25 (X-ray) 


a Ori 


19 


Full catalogue 






20 


Members 


11 (X-ray) 




21 


Members 


12 (Spectroscopy) 


IC348 


22 


Full catalogue 






23 


Members 


10 (Periodic Variables) 




24 


Members 


16 (X-ray) 




25 


Members 


15 (X-ray) 




26 


Members 


9 (Spectroscopy) 




27 


Members 


17 (Ha) 


h and \ P er 


28 


Full catalogue (combined) 






29 


Positionally isolated 





Table 2. The catalogues presented in this paper. 
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Figure 1. The difference in V — I for Cep OB3b between data 
taken in the system used here and using Cousins I filter. The line 
is the best fit through the data. 



V-I c - As demonstrated by the extreme example of CepOB3b 
(Av — 3 mags), even at a V-I — 3 the shift is < 0.1 mags; 
though note that the reddest standard in our V-I calibration 
is 2.87. This means over the colour range covered within this 
paper the effect is negligible. 



3 SEQUENCE ISOLATION 
AND FITTING 

For each field-of-view (FOV) the CMD includes a high pro- 
portion of background and foreground objects, therefore we 
faced the common problem of identifying members. For most 
fields we used literature sources for members. We discuss 
these fields in Sections l3.ll and l3.2l For h and x P er no mem- 
bership data were available for the group, or none at faint 
enough magnitudes. Here we can increase our confidence in 
selecting the correct sequence by selecting a limited area of 
the sky. This method will clearly include a higher proportion 
of non-members. The field treated in this way is discussed 
in Section 1331 



3.1 Sequences with memberships 

In Appendix [A] we collect literature membership data for 
each field, a summary of which is presented in Table [3] We 
use a wide range of membership criteria, each of which will 
have an implicit bias. As we are interested in the form of 
the sequences in CMD space it is important to be aware 
of and limit any bias incurred with respect to the CMD 
space. X-ray selection will be biased to weak lined T Tauri 
stars (WTTS) and binary systems l|Preibisch fc Feigelsonl 
120051 ). Variability exhibits a bias towards accreting objects, 
though periodic variability may exclude many classical T 
Tauri stars (CTTS) objects when observed over a long base- 
line. However periodic varia bility deduced from a short base- 
line does not carry this bias (jLittlefair et alj|2005f ). The final 
frequently used data source was Ha, which is generally bi- 
ased towards CTTS. 

There is also a more subtle form of bias due to pre- 
selection. Because in many cases surveying all objects in 



a field is not possible, an observer will necessarily select a 
subset of objects. Often this initial selection is based on the 
positions of objects in CMD space thus introducing a bias in 
CMD space. This is particularly applicable to spectroscopic 
data. Where this form of sample selection has occurred the 
selection criteria are stated clearly in Appendix |A1 

Samples derived from non-spectroscopic methods of 
membership selection have a higher probability of containing 
contamination from background (BG) or foreground (FG) 
objects. This leads to the need for colour-magnitude selec- 
tions to remove data points which lie far from the group 
sequence. Where we have done this the number of discarded 
objects is small relative to the total number of objects and 
has little effect on the final sequence (see Section 1 3 . 2 p . In 
each section where such a selection was used, it is shown in 
the CMD and its effect on the result is discussed. 



3.2 Fitting procedure 

To explain the stages involved in producing the empirical 
isochrones we use an example cluster, NGC2264. This clus- 
ter is a good illustration as it has a large number of members 
from various sources. 



3.2.1 NGC2264: an example 

Once the sequence was finalised to stars having properties 
indicative of cluster membership (see Appendix [XT} we fit- 
ted the median filtered sequence with a smooth curve. The 
median filter was used for two reasons. The median is neces- 
sarily tied to the values of real stars. In addition the median 
allows the effect of stars lying far from the sequence to be 
limited, as not all of our memberships are certain, for in- 
stance the X-ray sources in this group. The median of any 
sequence was calculated by binning the stars in V-band in- 
tervals and assigning the median V-I and V values to the 
median star for that bin. The bin sizes were tailored to each 
field, and we enforced a minimum number of stars in each 
bin. A cubic spline was fitted through the resulting list of 
median stars, over the colour range of the stars selected, with 
the gradients at the ends tied to a quadratic. The resulting 
fit for NGC2264 is shown in Figured 

To remove BG and FG objects and spurious matches 
which lie below the PMS a colour-magnitude selection was 
applied before fitting. We removed objects lying below the 
dotted line in Figure |3] We were concerned by the effect of 
this colour-magnitude selection on the spline functions, so 
we examined fits using data with and without a selection. 
The effect of removing the colour-magnitude selection for 
NGC2264 is shown in Figure El It has little effect on the 
resulting spline. 

The resulting fit is clearly an average through the spread 
in the colour mag nitude diagram (sometime s referred to as 
an age spread, see iBurningham et al. I l2005al , and references 
therein). Part of this spread will be due to unresolved bi- 
naries which lie above the single star sequence, lifting the 
spline above the single star sequence. Thus our splines do 
not necessarily follow the single-star sequence, although the 
displacement due to binary should be the same for all fields. 
Furthermore the scale of this effect is small, as can be seen 
in the CMDs were there is a clear separation between the 
binary and single star sequence, e.g NGC2547 (Figure H). 



Empirical isochrones 5 



Group 


Data Type 


Source 


Matching 


Offsets 








radius 


RA, DEC (arcsec) 


NOC2547 


X-rav ('ROSAT') 


I 


From 1 







Spectroscopy 


2 


1 arcsec 


o n 
u u 


NGC7160 


Extinction 


3 


1 arcsec 





NGC2264 


X-ray (ROSAT) 


4 


6 arcsec 







Ha 


5 


1 







Periodic Variables 


o 


1 arcsec 


u u 


Cep OB3b 


X ray (ROSAT) 


7 


From 7 otherwise 14 arcsec 


(both) 




"V raw ( HH A NT)R A APT^ 


9^ 




n n 
u u 






g 


1 iAvr'Uf^f I nnt li 1 

_1_ (XI Lo\_.L l l_J W Li 1 1 y 






Ha 


19 


1_ arcsec 









18 


_1_ (XI LSCL 





a Ori 


X-ray (XMM-NEWTON) 


11 


6 arcsec 







Spectroscopy 


12 


1 arcsec (both) 


(both) 


ONC 


Periodic Variables 


13 


1 arcsec 


+1.5 -0.25 and 




X-ray {CHANDRA) 


14 


2 arcsec 





IC348 


Periodic Variables 


10 


1 arcsec 


-0.1 -0.5 




X-ray {CHANDRA ACIS) 


16 


1 arcsec 


-0.5 




X-rav ("ROSAT") 


15 


14 arcsec 







Spectroscopy 


9 


1 arcsec (both) 


and +0.2 -0.5 




Ha 


17 


1 arcsec 


+0.2 -0.5 


NGC2362 


Ha (Li) 


20 


Members from 20 




A Ori 


Spectroscopy (Li) 


21 


Members from 21 




IC5146 


Ha 


22 


1 





NGC6530 


Ha 


23 


Members from 23 






X-ray {CHANDRA ACIS) 


24 


Members from 24 





Table 3. Table of lite rat ure sources and astrometr ic matching criteria. 1 [Jeffries fc Tollev l rtl998T), 2 
I Jeffries fe Oliveiral j2005h . 3 lSicilia-Aeuilar et al.l d2005Tl , 4 lFlaccomio et al.l lll999Tl. 5|Dahm fc Simonl d2005Tl . 
6 iLamm et al.l d2004fTj Secon d ROSAT PSPC catalo gue and iNaylor fc Fabianl lll999l). 8 iPozzo et atl 
J2003l) and lOgura et aLI (12002ft. 9 iLuhman et all (l2003h a nd iHerbi el lll998l). 10 ICohen et al. I d2004h and 
iLittlefair et al.l ^OOslT^l l ISanz-Forcada et al] ||2004|)~12 iBurningham et all l|2005bl l and iKenvon et al.l 



ll2005h. 13 iHerbst et al] 1120021 ). 14 iFlaccomio et al.l ||2003|) , 15 Second ROSAT PSPC c atalogue, 16 
IPreibisch fc Zinneckerl J2002I1 17 iHerbid Jl99STl. 18 Lit tlefair et al (in prepa ratio nl . 19 lOgura. et al] 
| |2002|). 20lDahml d2005T). 21 iDolan fc Mathieul J2OOII) . 22 iHerbig fc Dahml llioM ) . 23 ISung et al] (120001 1. 



(2002), 20 Dahm (2001 
24 |Prisinzano et al] l| 200 



21 IDolan fc Mathieul d200ll) , 22 IHerbig fc Dahml (120021) . 
25 ICetman et al] I^Ood Tand 26 lMendoza V. fc Gomel <|l98oT v 



Sometimes the spline fit does not precisely follow the 
sequence the eye picks out, for example 0.5 < V-I < 1 for 
NGC2264. Whether the spline or the eye is correct turns out 
to be immaterial. The deviations are never large enough to 
affect our conclusions. 



3.2.2 NGC2547 

The o bservations of NGC2547 were taken from lNavlor et all 
(2002), using Cousins / (/ c ). The photomet r y was calibrated 
with coefficients detailed in iNavlor et all $2002) with the 
reddest standard having a V-I of 2.7. The resulting fit, after 
sequence selection is shown in Figure [4] A colour-magnitude 
selection was used to clip out probable non-members. Only 
eight objects were removed, all lying far from the sequence 
and within the contamination. The fitted line lies slightly 
above the single star sequence, an effect due to binary stars. 



3.2.3 ONC 



The photometry (usin g I c ) and members hip probabilities for 
this region come from lffillenbrandl (|l997n . For the ONC the 
extinctions are calculated individually for each star, there- 
fore the order in which the fitting and conversion to absolute 
colours and magnitudes is carried out is crucial. Where av- 
erage extinctions have been used the entire space is simply 
translated and the fitted spline curve will not change shape. 
If however individual extinctions are used shifts may there- 
fore be different for each object, changing the shape of the se- 
quence. Consequently, for the ONC the extinctions must be 
applied before the fitting procedure. This is also the case for 
the later fields of IC348 and IC5146. The resulting sequence 
(extinction corrected) with additional membership informa- 
tion fitted is shown as Figure [5] All the stars shown were 
fitted to yield the spline curve shown. No colour-magnitude 
selection was utilised here as nearly all the stars within 15' 
of the centre are members of the ONC. This is because the 
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Figure 2. Stars selecte d as members of NG C2264. Circles are the 
periodic var iables from [ Lamm et al ] [|2004h . asterisks are X-ray 
sources fromlFlaccomio et a.1.1 Jl999T) . crosses are Ha sources from 
iDahm fe Simonl d200^" and t riang les are proper motion members 
from lMendoza V. fc Gomej dl980ri . The spline fit to the sequence 
after a colour-magnitude selection is shown. 




Figure 4. The full catalogue for the NGC2547 field (do ts). As- 
terisks are the X-ray sources fro m! Jeffries fc Tollevl dl998f) ■ circles 
are spectroscopic members from IjeffrieT'fe'cjliveiral (120031 . The 
X-ray sources, some of which clearly lie in the BG MS contamina- 
tion were removed using a colour-magnitude cut prior to fitting, 
shown here as the dotted blue line. 




Figure 3. The effect of choosing a wider colour-magnitude se- 
lection in NGC2264. The effect is limited by the large number of 
sources within the sequence. The solid curve is the result using no 
selection, with the dotted line being the colour-magnitude selec- 
tion and the dashed line the spline fit after the colour-magnitude 
selection is enforced i.e. only stars above the line were used to 
produce the spline curve. 



background stars are largely obscured by the dense cloud 
cloud behind the cluster. 



3.24 NGC7160 

The sequence and fit is displayed in Figure [6] The mem- 
bers, although sparse, represent a relatively unbiased sam- 
ple in CMD space due to the 'wide' criteria applied to assign 
memberships (see Appendix I A4[l , 



3.2.5 a On 

For a Ori a colour-magnitude selection was used to clip out 
the X-ray sources lying within the contamination. The re- 



in 



15 



25 





1 




I 












\ * 
































* 




* 


















£*: ° 






























o. 












1 




1 




i 








2 



Figure 5. The ONC photometry from lffillenbrandl 1 19971) (dots'). 
Circles are objects which are periodic variables from Hcrb st et al 



2002) . The asterisks are X-ray sources from iFlaccomio et al 

2003) . The best fit spline through all data points is also shown. 



suiting line fit is displayed in Figure [7] The applied colour- 
magnitude selection removed five objects all far from the 
sequence and within the contamination. 

3.2.6 Cep OB3b 

The selected members and fit for Cep OB3b are shown as 
Figure [8] A colour-magnitude selection has again been ap- 
plied, although as in Section 13.21 since most objects lie in 
the sequence the effect on the fit is negligible. 



3.2.7 IC348 

For IC348 individual extinctions were available for all but 
21 of the members. Therefore only those members with ex- 
tinctions have been used. A colour-magnitude selection has 
been applied removing five stars below the sequence. The 
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Figure 6. The full catalog ue for the NGC7160 field ( dots). 
Circles are members fro m ISicilia-Aguilar et al.l l|2004l ) and 
ISicilia-Aguilar et alj J2005h . 




Figu re 7. Stars selected as members in a Ori. Circles are 
iBurning ham et al.| |2005bl) memb ers. Asterisks are X-ray sources 
from ISanz-Forcada et alj (2001). Crosses are members from 
iKenvon et al.1 1|2005|) . The fitted spline curve is shown. The colour- 
magnitude selection is shown as the dotted line. 



individual dereddenings were applied prior to fitting. The 
dereddened sequence with the fitted curve can be seen as 
Figure [5] 




Figure 8. Stars selecte d as members in Cep OB3b. Asterisks 
are X -ray sources from iNavlor fc Fabianl jl999h . | Gctma n et alj 

and the second PSPC catalogue. Triangles are members 
fromlPozzo et alj<2003fl . Circles are Hct sources from Ogura ct al. 

Crosses are the periodic variables from Littlefair et al (in 
preparation). The fitted spline curve is shown. 




Figure 9. Stars se lected as members in IC348. T he asterisks are 
X-ray sources from lPreibisch Sz Zinneckerl (2002) and the Second 
ROSAT PSPC cata logue . Circles are the perio dic variables from 
ICohen et all J2004h and iLittlefair et alj J2005I) , Crosses are Ha 
sources from lHerbig] lll998h. Triangles are sp ectro scopic member s 
with extinctions from iLuhman et alj ||2003!) and [ Herbigl jl998l). 
Indivi dual extinctions from lLuhman~et**aL I (I2003T ) and iHerbiel 
(1998) have been applied before fitting. 



3.2.8 A On, NGC2362, IC5146 and NGC6530 

The data for A Ori are from lDolan fc Mathieul (120011 ) taken 
in the Johnson-Kron- Cousins system and calibrated to Lan- 
dolt standards, a colour-magnitude selection has been ap- 
plied removing three stars lying well a b ove t he sequence. 
The data for NGC2362 are from iDahml (|2005l ). taken with 
an I c filter and tied to Landolt standards. The Cousin s 
photometry for IC5146 comes from iHerbig fc Dahml (|2002l ). 
here only the stars shown as asterisks in Figure [12] were 
fitted. This contains photometry from two areas approx- 
imately 10' apar t . The photometry for NGC6530 is from 
iPrisinzano et all (|2005h . The members and spline fits for 
these fields can be seen as Figures [101 1111 [T2l and 1131 



3.3 Isolating the sequence for h and \ P er 

In this field there was no membership data of the form used 
above, so colour-magnitude selections were employed. In ad- 
dition as the group in question is clustered in the FOV it was 
possible to remove a proportion of FG and BG objects by 
selection of a limited area of the FOV. We however needed 
a guiding mechanism to allow us to pick the approximate 
central coordinates and a subsequent area to isolate. The 
following procedure was undertaken to isolate the cluster in 
celestial coordinates. 

(i) Lead Stars. We show our CMD for the full h and \ 
Per catalogue in Figure 1141 The sequence is clearly visible 
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> 15 




Figure 10. St ars selected as members in A Ori. The dots are Li 
members from lDolan fc Ma thicu (200l]). 




Figur e 13. The photometry of NGC6530 from [Prisinzanc^etjil 
12005) (dots). Crosses are X-ray sour ces from Prisinza no et al 
I 2005). Circles are weak Ha s tars fromlSung et al.l ||2000| ). Aster- 
isks are strong Ha stars from lSung et a.1.1 [|2000l V 




h and x Per 



Figure 11. Stars selected as members in NGC2362. The d ots are 
Ha an d spectroscopically confirmed PMS members from iDahml 
j2005h . 




Figure 12. Stars sel ected as members in IC5 146. Dots are likely 
PMS member s from [ Hcrbig & Dahm (2002). Asterisks are Ha 
members from|Herbie & Dahm ( 2002) with known spectral types 
excluding stars lying below the Pleiades MS. Only the asterisks 
were fitted. 




Figure 14. The full catalogue for the h and \ P er field (dots). 
The dotted line shows a colour-magnitude selection of V < 14th, 
showing the sequence clear of the contamination. The solid line 
is the best fitting spline to the selected sequence corrected to the 
mean extinction for the two clusters. 



and the brighter stars clearly lie blue-ward of the contam- 
ination. Therefore the positions of these stars on the sky 
were plotted, allowing us to trace the outline of the clusters. 
This can be seen in Figure IT5l 

Three other set s of lead star candidates were explored. 
lUribe et all 120021); proper m o tion m emb ers or photometri c 
members from ISlesnick et al.l (|2002l ); and iMuminovl i|l996l ). 
The coordinates from each set were carried forward through 
the following steps. However the central coordinates for all 
cases were very similar. 

(ii) Positional selection. A circular positional selection 
was applied to the full catalogue around the central coordi- 
nates, with an inclusion radius varying from 5' down to 1' 
in increments of 1'. The CMDs for each were examined and 
a 'by eye' selection made for the CMD yielding the clearest 
sequence, optimising the member to field-object ratio. The 
best results were obtained with a radius of 2'. The final cen- 
tral coordinates for each cluster are: a— 2 22' 5.02" S=+57 
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2.43 2.40 2.38 2.36 2.34 2.32 2.30 2.28 
Bight Ascension 

Figure 15. The positions of stars in the h and \ Per catalogue. 
Asterisks are stars with V < 14th, and the large circles show the 
regions we selected. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

V-I 



Figure 16. Stars within 2 arcmin radius of the centres of h and 
X Per. Stars in the h Per region have been shifted in extinction 
to match those \ P er - The dotted lines are the colour-magnitude 
cuts used, and the solid curve the best fitting spline. 



T 43.44" (x Per), q=2 18' 58.76" <5=±57 8' 16.54" (h Per), 
J2000. 

Since the measured ages in virtually all the literature for 
h and x P er are the same, it was possible to increase our 
confidence of selecting the correct sequence in CMD space 
by combining the CMDs from both clusters. To achieve this 
any differences in extinction and distance modulus must be 
accounted for. We derived extinctions and distance modulus 
shifts in the sense h Per minus x Per using three methods. 

(a) Matching the curve in the MS at V=16 mags for both 
clusters by eye, we found A v = 0.08±0.05, d m = 0.07±0.05 
mags. 

(b) Matching the centre of the UMS (Upper Main Sequence) 
of x Per to the centre of the more scattered UMS in h Per 
gave a shift in \d m \ of < 0.05, with a shift in extinction of 
A v =0.2 ±0.05. 

(c) The Q-method (see Section [47T]) gave A v = 0.157 ± 
0.0212. 

We chose the Q-method value for the extinction and no shift 
in dm ■ The scatter in the different derivations suggests an 
uncertainty of approximately ±0.05 in extinction and dis- 
tance modulus. This shift was applied to h Per, and the 
combined catalogue used from thereon (Table 28) . The final 
area selected can be seen overlaid on Figure [TBI 

(iii) Colour-magnitude selection and fitting. The fi- 
nal sequence was then fitted as described in Section 13.21 
Cuts in colour-magnitude space were then used to clip out 
stars lying far from the sequence. The cuts were varied until 
a best fit of the spline curve to the sequence was achieved. 
The final sequence, cuts and spline curve are shown in Fig- 
ure [TB] To further increase our confidence in selection of the 
correct sequence the best fit spline was shifted to the mean 
extinction for both clusters and overlaid on the CMD show- 
ing the full catalogue. The spline curve clearly follows the 
sequence one would select by eye from this CMD, see Figure 

M 



4 PARAMETERS 

To compare the data for our sequences on a single CMD we 
required estimates for the extinction (Ay) and true distance 
modulus (d m ). To obtain a self-consistent set of empirical 
isochrones it is clearly better to use a single method to de- 
termine extinction. This was possible for a range of fields 
using the Q-method as detailed below. Where data for the 
Q-method were unavailable we have used values from the 
literature, given in Appendix [B] A summary of the distance 
moduli and extinctions we used can be found in Table IB2I 

4.1 Q-method 

Using UBV photometry from the litera ture we applied the 
Q-method of IJohnson fc Morg an (1953). This method uses 
the colour relationships of OB s tars and is only val i d over 
a small range of spectral types. IJohnson fc Morgan! (|l953T ) 
quote (B- V) = -0.009±0.337Q to be valid for -0.80 < Q < 
—0.05, converting to an intrinsic colour range of —0.28 <(B- 
V)o< —0.03. This was used in conjunction with the initial 
parameter estimates to calculate an apparent colour range 
over which the method was applicable in each field. We then 
derived a mean extinction which we could in turn apply to 
our selection. The loop of applying a mean extinction, cal- 
culating Q from the stars within the colour range and calcu- 
lating a new mean extinction was iterated until the number 
of stars used in successive iterations remained constant. The 
RMS around the mean value divided by the square root of 
the number of points provides an uncertainty. 

4-1.1 h and x P er 

Here UB V photometry was obtained from ISlesnick et alJ 
(2002). In Section 13.31 we limited our study to two circu- 
lar fields centred on the clusters. As extinction varies across 
the field of view for these two clus t ers, w e have limited our- 
selves to stars from lSlesnick et alJ |2002l ) lying within these 
areas. The result is shown as Figure [T71 The difference of the 
median values of extinction using Q between the two posi- 
tionally selected areas was found. The shift derived for h 
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Figure 17. Le/i. The sequences for h and \ P er O and B stars 
after the application of the Q-mcthod. Right. The extinction. The 
circles are \ Per and the crosses are h Per. The resulting median 
values of Ay=l.b71 and 1.728 for \ an d h Per respectively are 
shown. 



Per was Av = 0.157 which has already been applied in Sec- 
tion [3]3] Therefore the reddening to \ Per only was needed, 
being A v = 1.571 ± 0.079. 



4.1.2 NGC2547, NG2264 and Cep OB 3b 

For NGC2547 we used the photometry of IClarial (|l982f ). 
The resulting median extinction value was 0.157 ± 0.008 
mags. For NGC2264 UBV member photometry from 
IMendoza V. fc Gomez! (|l980h yielded a median Av = 
0.371 ± 0.0409. Th e final field i s Cep OB3b, where the pho- 
tometry was from Ijordi et al] l|l99fj| ). with a resulting me- 
dian value of A v = 2.88 ± 0.039. 



5 DATA COMPARISON 

These data present an opportunity to make a qualitative 
assessment of how the PMS changes with age. To enable 
us to do this the sequences were shifted by the adopted 
distance moduli and extinctions (see Figure [T8)l. We selected 
four of our best sequences spanning as large an age range 
as possible. Data for NGC2264 has been supp lemented by 
photometry from IMendoza V. fc Gomez! §980). They used 
radial velocity to assign membership probabilities, and we 
plot only stars with membership probabilities of > 90%. 
This allowed us to plot the sequence members above Vb=ll, 
where the stars are saturated in our catalogue. 



5.1 Sequence spreads 

The sequences in Figure [TH] appear to "crystallise" with age. 
The absolute spread in CMD space of the sequences reduces 
as the age of the field increases. The trend appears to con- 
tinue in the fields older than those presented in Figure [TBI 
with NGC2547 (~30 Myrs) showing a particularly clear se- 
quence. Sequences older than 5 Myrs do not appear to show 
large spreads. However, sequences younger than this show 
large spreads not explicable by photometric errors, variabil- 
ity and binarity, as shown in the case of Cep OB3b and a Ori 



bv lBurningham et al] l|2005ah . This contrast with the Upper 
Sco-Cen OB association at ~ 5—6 Myrs, whi ch has an appar- 
ently wholly explicable photometric spread dPreibisch et al] 
2002). It appears there is a change at about 5 Myrs. 

This suggests a possible cause of the anomalous spreads; 
namely accretion or infall history. The half-life for accretion 
discs is estimated to be 3 Myrs, therefore for low-mass ob- 
jects o ne would expect the discs to have disapp eared by ~ 6 
Myrs (|Qliveira et al] 120061 ; Irlaisch et all 120011 ). Accretion, 
from a disc increases the continuum luminosity along with 
various other effects, possibly pushing the stars blue-ward, 
resulting in a sc atter both temporally older and younger 
l|Tout et al]|l999h . This component of the spread would not 
be present in sequences older than ~ 6 Myrs. It should be 
noted however that this half-life is derived from the presence 
of dust; the half-life of gas within these discs is less certain. 

An obvious alternative is that the groups are not co- 
eval. The timescale of star formation itself must produce an 
intrinsic spread in age and hence in position in the CMD. In 
addition there is strong evidence for the episodic star forma- 
tion within a group, sub-group, cluster or association. Thus 
another explanation for the secular evolution of the sequence 
spreads we observe is in Figure [18] could be an initial age 
spread evolving to older ages and the absolute spread there- 
fore reducing as the isochrones bunch up in CMD space, 
i.e. d(age)/dV increases. If the spread within the CMD is 
due to a genuine age spread, since our empirical isochrones 
are based on median magnitudes they will be biased toward 
younger ages. 



5.2 Evidence for a gap in a V, V-I CMD 

The ONC has a discernible gap between the PMS stars 
and the MS stars, at the point where the PMS isochrone 
meets the zero age main sequences (ZAMS). This can be 
seen clearly in Figure 1181 There also appears to be a gap, 
again where the PMS meets the MS in h and x P er an d 
NGC2264 (Figure fig)) . A gap is evident in the photometry 
for NGC6530 and perhaps IC5146 (see Figures [T3l and [T2")) . 
A further excellent e xample of the gap can be seen in the 
CMD of NGC4755 in lLvra et ail l|2006h . 

To ascertain whether this gap between MS and PMS 
stars is a general feature we used the isochrones to find 
the position (in CMD space) of the PMS-MS connection. 
In some of our younger fields the gap is situated around our 
bright magnitude limit. In these cases we cannot be sure of 
the gap's existence. There is however a discernible dearth of 
members at the head of the PMS, which can be seen in Cep 
OB3b, a Ori, IC348, NGC2362 and A Ori. This explains the 
general shape of many CMDs of PMS regions, which appear 
to consist of a sequence parallel and redward of the con- 
tamination which disappears at brighter magnitudes. The 
brightest stars at the head of this sequence must represent 
the top of the PMS, with the MS appearing blue- ward of 
this. The MS, and therefore the gap, is often obscured by 
contamination. 

We could not definitively detect the gap in NGC7160 
or NGC2547. In the case of NGC7160 the gap would fall in 
an area devoid of membership information. Finally for our 
oldest field NGC2547 the PMS and MS isochrones lie very 
close together meaning the gap, if present, would be very 
small compared to the younger fields. 
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It is apparent that this phenomena occurs in all datasets 
where we would expect to observe it. Th is feature del i neates 
the transition from PMS to MS, hence IStolte et all (|2004l ) 
call this the PMS/MS transition region (they identified it 
in a J 3 , J s - K a CMD of NGC3603). We prefer the term 
radiative convective (R-C) gap as we feel this emphasizes 
the physics involved which we outline in the next section. 



5.3 The origin of R-C gap 

From the theoretical isochrones an explanation for this gap 
becomes apparent (see Figure [19}. Physically the R-C gap is 
the transition phase of stars from a fully convective object 
to the development of a radiative core. The Hayashi tracks 
which the lower mass PMS stars descend in CMD space 
represent the gravitational contraction of a fully convective 
object. As the star heats and the core density gradually 
increases, at some stage the core becomes radiative, and the 
track changes to a Henyey track which is almost horizontal 
in a V/V-I CMD. For massive stars this h appens earlier in 
the star's evolution than for low-mass stars (jCollinsll 19891 ); a 
7Mq star will develop a radiative core almost immediately. 

The development and evolution of the radiative core 
causes the separation of the PMS and MS. As the transi- 
tion proceeds the mass tracks are almost parallel with the 
isochrones, meaning evolution across the area of CMD space 
between the PMS and MS is rapid compared to the descent 
down the Hayashi tracks. The difference in motion in CMD 
space is clearly demonstrated in Figure 1191 Here the dis- 
tance between a 1 Myr and 3 Myr isochrone for a star on a 
fully convective track is small when compared to a star at 
the head of the PMS which will have joined the MS by 3 
Myrs, moving much farther in CMD space. This means at 
any given age a sequence should have a relative sparsity of 
stars in this gap compared to the slowly evolving PMS and 
almost stationary MS (on these timescales). 

As is visible from the data in Figure [TH] and the theory 
in Figure [19] the size of the dislocation between the MS and 
PMS is a function of age. For the younger fields the PMS lies 
far from the MS and the transition occurs at high masses. 
The higher the mass of a star the earlier the formation of a 
radiative core and the sooner it changes to a Henyey track. 
Higher mass objects move to the MS more swiftly along 
these tracks. The reverse is also true. As the field becomes 
older, the mass at which the core develops falls and the 
evolution of the star progresses more slowly. 

The age dependency of the size of the gap has impor- 
tant ramifications. If it is possible to identify the head of the 
PMS and the MS members for a given field, the separation 
of these will be a distance and reddening independent mea- 
sure of the age. There are several restrictions. If the field 
is older than around 15 Myrs the gap becomes too small 
to detect reliably. The gap may also lie within an area of 
high contamination, meaning many objects will need to be 
assessed for membership to be certain of detecting the gap. 
In some cases variable extinction will prevent photometric 
detection of the gap, as the relative positions of the stars 
shift in a CMD; in these cases individual extinctions will 
need to be found. As the size of the gap depends critically 
on the formation and growth of the radiative core, it may 
also represent an excellent test of the stellar interior models. 



T 
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Figure 19. Isochrones (solid blue lines) from Siess ct al. I teooch 
for the approximate ages of the four fields in Figure [T8l 1, 3 and 
13 Myrs. Mass tracks (dashed lines) are show for 7, 6, 5, 4, 3, 2, 
1.2, 1.0, and O.SMq. The evolution between 1 and 3 Myrs of 3 and 
IMq stars are shown as filled circles, with the track highlighted 
in red. 

6 RELATIVE AGES 

The sequences displayed in Figure [18] are the results af- 
ter membership selection. To create a relative age ladder 
however the sequences must appear on the same CMD, su- 
perimposed over one another. Such a figure would be con- 
fused due to the number of points and the spreads of the 
sequences. Therefore at this stage we move to representing 
the sequences by the spline fits. To aid analysis we also cre- 
at ed a ZAMS subtr acted CMD, using the ZAMS relation 
of Siess et all (2000). In these diagrams the points on the 
spline fits each have the V-I of the MS star at equivalent 
V subtracted from their V-I. This separates the sequences 
much more clearly. Figure[20]shows a sample of the empirical 
isochrones. 

The disadvantage of the ZAMS subtraction process is 
that it enhances the residuals of the spline fitting procedure. 
As can be seen in Figure[5T]at the R-C gap there is a dearth 
of data which causes the spline fit to degrade. This also 
causes the fit to the MS after the R-C gap to be poor as 
the spline fits enforce a gradient continuity criterion. These 
residuals have been acceptable (if noted) until now but are 
greatly magnified in the ZAMS subtracted space. Therefore 
we have cut the splines in ZAMS subtracted space at a Vo = 
3. 

6.1 Theory comparison 

It is useful at this stage to compare the empirical isochrones 
we have created to their theo retical counterpart s. Although 
we have used isochrones from Siess et all (|200Ch , the effects 
described are apparent with all the other isochrones we have 
examined, though a complete investigation of all the models 
lies outside the scope of this paper. 

We have adopted an ap proximately solar metallicity, 
Z = 0.02. Ijames et alj ([2006) shows that for a a number of 
star forming regions the metallicty is slightly sub-solar, it is 
also likely that different clusters, associations, groups or sub- 
groups will have differing compositions. Therefore it is perti- 
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Figure 18. The likely PMS objects for fields of different ages, plotted in dereddened colours and absolute magnitudes. The spread in 
absolute co lour magnitude space clearly decreases as age increases. The ZAMS and isochrones of approximately the correct age from 
ISiess et al. I (l2000f) are shown. 



nent to ask what age difference would be obtained by fitting 
a higher or lower metallicity isochrone to our representa- 
tive sequence of 3 Myrs. If we compare a 3 Myr Z = 0.02 
isochrone with isochrones of Z — 0.01 and Z — 0.04 the 
closest matches have approximate ages of 2.2 and 3.5 Myrs 
respectively. Thus changes in the metallicity by a factor two 
are too small to affect the discussion which follows. 

In Figure [20l we can see a minor problem. NGC2547 can 
be seen to cut across the younger sequences, this is probably 
due to observations being made using a different / filter (see 
Section [2]3} coupled with the fact that the reddest standard 
is at a V-I of 2.5 in this case, exaggerated by the ZAMS 
subtraction process. 

In addition, our comparison reveals three serious con- 
cerns, (i) Figures [181 and [2l1 show that whilst an acceptable 
fit to the theoretical isochrones can be obtained over most of 
the sequence, the model isochrone lies above the data at the 
boundaries of the R-C gap for NGC2264 and h and \ P er - 
Furthermore the MS members overlap the PMS in magni- 
tude space for all our fields, in direct contradiction to the 
isochrones, assuming a coeval population (see Section[9]for a 
further discussion of this) . (ii) For younger sequences a the- 
oretical isochrone fitted to the PMS is a poor match to the 
MS. (iii) Although the young sequences have a large scat- 
ter, the PMS isochrones fail to consistently follow the centre 




Figure 20. A sample of the empirical isochrones. Youngest to 
oldest, the ONC, NGC7160, h and \ Per and NGC2547. The left 
panel is the empirical isochrones in an abso lute V, V-I CMP , 
with the lowest curve the ZAMS relation of Siess et al ] l|200d) . 
The right panel is the ZAMS subtracted space; at a given V the 
ZAMS V-I is subtracted from the V-I of the spline. The vertical 
dashed line shows the position of the ZAMS 




of this spread. For example, in the NGC2264 most stars lie 
above the isochrone at V-I> 1, but below it at V-I< 1. 
The conclusion is that when fitting isochrones to sequences 
across a limited colour range, one must be aware that the fit 
will not necessarily be good across the rest of the sequence. 
This implies that the age derived depends on which part of 
the sequence is fitted. This is why empirical isochrones are 
potentially superior to their theoretical counterparts. 

6.2 Selecting the fiducial sequences 

As we shall see later, we can roughly group the youngest 
objects into three groups with ages of 1, 3 and 5 Myrs. 
The sequences older than this are much more clearly sep- 
arated. These three groups are separated in age such that 
using the maximum uncertainty in distance modulus (typ- 
ically 0.2 mags) could shift a sequence by one group, older 
or younger. For use as fiducial sequences we select the spline 
curves of one sequence from each group, typifying the group. 
We preferentially select the sequence from each group with 
the most certain distance and age. These fiducial isochrones 
can then be used to derive relative ages to further fields, with 
selection of the fiducials guided by the sequence's placement 
in colour-magnitude space. The group at around 5 Myrs is 
problematic, as discussed later it contains only two clusters 
a Ori and 1C348; a Ori has an uncertain distance. However, 
IC348 is found to be much older (in this work) than the 
age given in the literature. Therefore a Ori has been cho- 
sen as the fiducial for this group, despite its large distance 
uncertainty. 

The fiducial empirical isochrones selected were the 
ONC, NGC2264, a Ori and h and \ Per, with a s sumed 
ages (from the literature or fitting of ISiess et"al] (|2000l ) 
isochrones) of 1, 3, 4 and 13 Myrs. These fiducials will be 
used in Section [7] The fiducials are shown in Figure [22] 



7 USING A RELATIVE AGE LADDER 

We can now use our age ladder to determine relative ages for 
several fields. Throughout the figures in this section the fidu- 



cial isochrones are marked as dashed lines, with the subject 
sequence appearing as a solid line. 



7.1 Cep OB3b 

For Cep OB3b the membership list has been discussed in 
Section IA6I and the CMD with the members overlaid is 
shown as Figure IA3I Estimates for the distance modulus 
are from Section TB7I and the extinction from Section [4.1.2l 
We use the spline fit to the data (Section 1 3 . 2 p , over- 
laid on the fiducial isochrones from Section 16.21 to age the 
subgroup. This reveals an interesting result. The subgroup 
Cep OB3b can clearly be seen to lie on the same locus as 
NGC2264, with a Ori being slightly older, implying an age 
for Cep OB3b of ~ 3 Myrs (Figure 1231. The b e st pr evi- 
ous age estimate for this subgroup is IJordi et al.l (|l996l ) at 
5.5 Myrs. The older estimate would have made Cep OB3b 
the oldest group to contain molecular material. We therefore 
conclude that the lower age of 3 Myrs is more likely. The dis- 
tance modulus uncertainty range is ±0.20 mags. The effect 
of this can be seen in Figure l2"4l where the relative position of 
Cep OB3b does not change dramatically even when the en- 
tire uncertainty budget is used. In addition the photometry 
shown in Figure [A"3l reveals evidence for an R-C gap at the 
head of the PMS in a position similar to that of NGC2264 
at an age of ~3 Myrs, also suggesting a similar age for the 
two groups. 



7.2 IC348 

Using our relative age system, as in Section 17.11 we found 
1C348 to have an age between that of a Ori and NGC2264. 
The empirical isochrone therefore suggests that the cluster 
is somewhat older than NGC2264 and indeed Cep OB3b. 
This can be seen in Figure 1251 This contrasts w ith the lit- 
erature age of 2-3 Myrs (e.g. IHaisch et al.ll20011 ). To obtain 
an age of 3 Myrs or less would require a shift of ~ 0.5 mags 
in the distance modulus, yet the uncertainties in the dis- 
tance modulus to IC348 are significantly smaller than this, 
-0.16/ + 0.14 mag (See Appendix [Bj . 
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Figure 21. The selected members for the ONC, Cep OB3b, h and \ Per and NGC2547. The solid line is the empirical isochronc. The 
dashed line is the best fit isochrone (by eye) from lSiess et all (2000) (~ 1.5, ~ 3, ~ 13 and ~ 16 Myrs respectively). The sequences and 
empirical isochrones can be seen to move away from the theoretical counterparts to varying degrees. 





1.0 1.2 1, 



Figure 24. The sequence for Cep OB3b assuming thee different 
distance moduli (solid lines). NGC2264, h and x Per, a Ori, and 
the ONC are shown as dashed lines for reference. 



Figure 25. The fitted sequence for IC348 (solid line), h and \ 
Per, a Ori, NGC2264 and the ONC are shown as dashed lines for 
reference. 



7.3 A Ori 

Figure l2"ul shows the fiducial sequence plot with A Ori placed 
at distance moduli of 7.77, 7.90 and 8.07 mags, represent- 
ing the limits of the HIPPARCOS measurement. Figure [26] 



shows A Ori is largely degenerate with NGC2264, and there- 
fore has an age of ~ 3 Myrs. This is certainly consistent with 
the MS turn-off age of 4 Myr de rived from UBV photom- 
etry by iMurdin fc Penstonl (|l977t ) but not with the age of 
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1.0 1.5 2.0 2.5 0.2 0.4 0.6 0.8 1.0 1.2 i. 

(V I), 

Figure 26. The sequence for A Ori assuming the dis t ance mod- 
uli from lDolan fc Mathieul l|200ll l. lMurdin fc Periston! ljl977h and 
HIPPARCOS (solid lines). The ONC, NGC2264, a Ori and h and 
X Per are shown as dashed lines for reference. 
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Figure 27. The fitted sequence for NGC2362 (solid line). The 
ONC, NGC2264, a Ori and h and \ P er are shown as dashed 
lines for reference. 



6-7 Myr derived by iDolan fc Mathieul (|200ll ) using narrow 
band photometry of high mass stars). 



7.4 NGC2362 

The relative age plot can be seen in Figure 1271 Here 
NGC2362 appears between NGC2264 and a Ori, giving it 
an age o f 3-4 Myrs. T his is consistent with the ages de- 
ri yed by iDahml (feoOrJ) of Myrs using the isochrones 

of lD'Antonafe MazziteUl dl997h . and 3.5-5 Myrs from the 
isochrones of [Baraffe et al.l (| 1998T ) . However, the shift in dis- 
tance modulus to achieve an age of greater than 4 Myrs 
would be ~ 0.4 mags; larger than the typical distance mod- 
ulus uncertainty. 



7.5 IC5146 

The age plot is displayed as Figure 1281 IC5146 appears to 
have an age approximately equal to that of the ONC. The 
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Figure 28. The fitted seque nce for IC5146 (solid lin e), assuming 
the distance modulus from iHerbig fc Dahml (|2002f l. The ONC, 
NGC2264, a Ori and h and \ P er are shown as dashed lines 
for reference. The vertical dashed line shows the position of the 
ZAMS 



median age found by IHerbig fc Dahml l)2002l ) is ~ 1 Myr, 
agreeing with our result. 



7.6 NGC6530 



Photo metry from this cluster is taken from lPrisinzano et al.l 
(2005). They only pro vide pho t omet ry of stars which are 
also in the catalogues of lWalkerl (Il957t) orlKilambil ll 19771 ) or 
matched to an X-ray source from Damiani et al. (2004), and 
are brighter than V = 1 8. IPrisinzano et al.l l|2005l) use the 
isochrones of lSiess et al.l (|2000|) fitted t o PMS stars to obtain 
a median age of 2.3 Myrs. Sung et all l|2000! ) find an age of 
1.5 Myrs from fitting isochrones to an Ha selected PM S after 
conversion to an H-R diagram. iDamiani et al.l II2004| ), using 
X-ray selection and the photometry of ISung et aiT (120001 ) 
find an age of 0.5-1.5 Myrs. As for previous fields this spline 
has been overlaid on the fiducial diagram s een as Figur e |2"91 
Two d istance moduli are displayed, that of IPrisinzano et all 
l|2005l ) (d m ~ 10.48 mags) and that of ISung et all (|2000h 
(dm — 11.25±0.1 mags). The sequence in Figure l2"91appears 
to lie above the ONC for the greater distance modulus and 
coincident, or perhaps just below, the ONC, for the lower 
distance modulus. Suggesting an age of < 1 Myr at the 
further distance or ~ 1 — 2 Myrs for the closer distance. 
Thus we cannot differentiate between the literature ages. 



7.7 NGC7160 and NGC2547 

These clusters have older age estimates than the fields aged 
using our earlier fiducial selection (Section 16 . 2f) . NGC7160 
has an age estimate of ~ 10 M yrs from isochrone fitting 
from ISicilia-Aguilar et all l|2005h . NGC2547 has age esti- 
m ates varying f rom 3 0-45 Myrs from MS isochrone fitting 
in lNav lor et al. ( 2002]) and the Lithium Depletion Boundary 
in I Jeffries fc Oliveiral (|2005l ). Therefore we plot these empir- 
ical isochrones with h and \ P er as our on ly useful fiducial. 
NGC7160 is clearly younger than h and % Per, making the 
current age estimates reasonable. NGC2547 cuts across the 
other sequences at redder V-I. This could be a problem with 
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Figure 29. The fitted sequence for NGC6530 (solid line). The 
ONC, NGC2264, a Ori and h and \ P er are shown as dashed 
lines for reference. 




Figure 30. The fitted sequences (solid lines) for (younger to 
older); NGC7160 (blue, ~ 10) and NGC2547 (red,~ 30-45 Myrs) 
. The dashed line is the spline curve for h and x P er (— 13 Myrs). 
The vertical dashed line shows the position of the ZAMS 



the standard calibrations, which extends only to a V-I of 2.5 
in the NGC2547 reduction. 



8 IMPLICATIONS OF THE REVISED AGES 

The consequences of our revision of the relative ages of 
IC348, NGC2264, NGC2362 and the ONC are particularly 
pertinent to recent studies in two areas. Firstly in gy- 
rochronology. Large-sample rotation studies now exist for a 
subset of young clusters and from these a conventional view 
of the rotational evolution of stars between 1 and 5 Myr is 
developing. The rotational period distribution of the ONC 
is strongly bi-modal. This distribution is believed to evolve 
into a uni-modal distribution such as that of NGC2264 and 
NGC2362. IC348 has a distribution of rotation rates simi- 
lar to that of the ONC, and yet we find it to be older than 
NGC2264. This clearly raises a concern for the universality 
of any evolutionary model. 

The second area our findings affect is the field of disc 
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Figure 31. Fraction of JHKL excess sources, used a disc indicator 
for th e ONC, IC348, NGC2264 and NGC2362 from lHaisch et al] 
(2001); left panel, against nominal a ges found in this work, 
right panel, against cluster ages used in H aisch et al ] ||200lft . The 
straight line is a simple linear regression fit. 



evolution. Irlaisch et all <|200lh plot the fraction of JHKL 
excess sources, used to infer a disc fraction as a function of 
age for different clusters or groups. This plot reveals a clear 
decrease in the inferred disc fraction with time. However, 
from this work two clusters, IC348 and NGC2264 should be 
relocated within this plot leading to a change in the finer 
structure of the trend at younger ages. The disc fractions 
for the ONC, IC348, NGC2264 and NGC2362 are ~ 80%, 
~ 65%, ~ 52% and ~ 12% respectively. The relative age 
order of these four clusters found in this work (youngest to 
oldest) is the ONC, NGC2264, NGC2362 and IC348. There- 
fore they do not show a secular decline in disc fraction. This 
is illustrated most clearly by the Figure |3"T1 showing the disc 
fraction as a function age from lHaisch et all (|200ll ) and as 
a function of the nominal ages found in this work. 

Interestingly, IC348 is out of place in rotation rate and 
disc evolution. If disc evolution drives the rotational evo- 
lution, say through disc locking, this is perhaps unsurpris- 
ing. The key here may be the local environment. There are 
no O-stars associated with IC348. NGC2362 however has 
6 O-stars stars designate d as cluster members, the ONC 
has 16 and NGC2264 13 (|Mai'z-Apellaniz et aQBooj ). This 
supports the idea that a larger fraction of the PMS stars 
in IC348 have retained discs due to the absence of winds 
and/or ionising radiation from massive stars. Our relative 
ages may, therefore, be amongst the first indications of the 
importance o f environment on d isc and rotation rate evolu- 
tion (see also lStolte et alj|2004l ). 

An a l ternat ive explanation has been put forward in 



Lyo et al 



Lvo et al 



2003) by indentifying NGC2362 as anomolous. 
(120031 ) suggest this could be due to an incorrect 



age determination; this is ruled out is by this work. They 
also suggest it could be an artifact of the poor limiting mag- 
nitude achieved when using the L band, meaning this band 
is very sensitive to distance and provides a disc fraction de- 
rived from a different mass range for clusters at an unequal 
distance. Therefore any conclusions reached on the secular 
trends of disc fractions must be tentative, due to the het- 
erogenous nature of the stellar mass range used. 
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9 CONCLUSIONS 

We have presented photometry for a large selection of young 
stellar populations in differing environments. The selec- 
tion criteria and some of the photometry is heterogeneous 
and from many literature sources. Despite these difficulties 
this attempt at producing a self-consistent set of empiri- 
cal isochrones has been largely successful. We have charac- 
terised and fitted the sequences in CMD space, and grouped 
the sequences by position in the CMD - effectively age. Then 
using the clearest sequences typifying each group, a relative 
age ladder was developed for use with additional sequences. 
Our final groupings for all the sequences we have data for 
were as follows. 

(i) ~1 Myr; The ONC, NGC6530 and !C5f46. 

(ii) ~3 Myrs; Cep OB3b, NGC2362, A Ori and NGC2264. 
(hi) ~4-5 Myrs; a Ori and 1C348. 

(iv) ~10 Myrs; NGC7f 60-uncertain due to lack of mem- 
bers. 

(v) ~13 Myrs; h and x Per. 

(vi) ~>30 Myrs; NGC2547 

The distance moduli remain the largest source of un- 
certainty, allowing the possibility of changing a sequence by 
plus or minus one group in age. Cep OB3b is thus some- 
what younger than previously thought, and 1C348 somewhat 
older. Both the disc fraction and rotation-period distribution 
for IC348 imply it is younger than our derived age. Since, in 
contrast to the other regions we have studied IC348 has no 
O-stars, this suggests that the winds and/or ionising radi- 
ation from these stars have a significant effect of accretion 
disc evolution. 

The data also allow us to follow the feature in the CMD 
we term the R-C gap. In most CMDs this feature appears as 
a gap between the PMS and MS members, or as a dearth of 
members at the head of the PMS in some younger sequences. 
We have been able to trace the evolution of the gap from 1 
Myr through to the intermediate age (in the context of these 
data) clusters of h and x Per. We have explored the physics 
behind the appearance of a gap in the V, V-I sequences and 
the evolution of the gap with age. It is encouraging that 
the size of the gap, where it is visible, follows the same age 
ordering as found using the empirical isochrones. Further to 
this we have touched upon the possibility of using the size 
of the R-C gap as a distance and extinction independent age 
indicator. 

The R-C gap could also be used as a powerful tool to 
calculate age spreads within a sequence. If one could ac- 
count for all possible forms of spread within a CMD, and 
had memberships and spectral types/masses for all stars at 
the boundaries of the R-C gap. By finding the lowest mass 
star to have already traversed the R-C gap, and then finding 
the youngest star (via isochrones) of the same mass on the 
PMS one could calculate an age spread. 

Within this study is a limited investigation of the 
fit across sequences of theoretical isochrones. We conclude 
that isochrone formulations are self consistent when applied 
within a restricted mass or colour range not bridging the 
R-C gap. They tend to become less reliable when we have 
attempted to fit both the MS and PMS i.e. both sides of 
the R-C gap, excluding the intermediate age clusters, h and 
X Per. Most recent studies in the literature concentrate on 



either the turn-on region of a sequence or the low-mass sec- 
tions, meaning that the inconsistency of the fits is not always 
apparent. 

Finally, we tentatively suggest a critical change in the 
nature of spreads within a CMD at ~5 Myrs. Above ~5 
Myrs the sequence members are generally not associated 
with molecular material and have wholly explicable photo- 
metric spreads. Below ~5 Myrs the sequence members are 
often associated with molecular material and show inexpli- 
cably large photometric spreads. 

To revisit these or other clusters with a larger photo- 
metric study applying the same mechanisms as in this paper 
would be beneficial and we believe could provide a bench- 
mark collection of empirical isochrones for comparison to 
theory. In unison with this a comprehensive study of the ap- 
plication and limitations of the R-C gap as an age indicator 
would be a logical next step. 
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A3 ONC 




Figure Al. The full catalogue for the NGC2264 field ( dots). Cir- 
cles are the periodic var iables from [ Lamm et al asterisks 
are X-ray sourc es from |Flaccom io et al, I jl99Sh and crosses are 
Ha sources from lDahm fc Simo nl (|2005h . 



APPENDIX A: LITERATURE MEMBERSHIPS 

This section details the memberships for each field. Table [3] 
shows the group, the data type and the reference for each of 
the groups in this section. 



iHillenbrandl (|l997l ) calculated extinctions for a subset of 
stars, using spectroscopic measurements of stars already 
identified as members from a range of literature sources, 
in a spatially selected region encompassing the inner 15' ra- 
dius. We have used only those stars with spectroscopically 
calculated extinctions with assigned members hip probabil- 
ities a bove 80%. Periodic vari ables were from Herbst et all 
(|2002h and X-ray sources from lFlaccomio et al.l (120031 ). 



A4 NGC7160 



Members were taken from ISicilia-Aguilar et all l|2004h and 
ISicilia-Aguilar et al.l l|2005h (Table 7) . In these studies a pho- 
tometric sample of stars was selected, initially by considering 
those lying above the 100 Myr isochrone in a V, V-I CMD 
and lying within a magnitude range of F=15-19 mags. A 
further selection was then made of stars with RI variability, 
then using [/-band and 2MASS colors (biased to accreting 
stars). Spectroscopy was then used to measure extinctions 
and therefore assign memberships. The probability of mem- 
bership was fixed to the distance in standard deviations (a) 
from the average extinction value for the group. This relies 
on intrinsic colours. The members selected are those within 
la of the average extinction. 



Al NGC2264 

The X-ray sources (Table 10) picked out a reasonably clear 
sequence (see Figure IA1|) . although there were also some 
X-ray active foreground and background ob jects. The spec- 
trosco pic members (Table 11) came from iDahm fc Simonl 
l|2005h who selected members using Ha strength from nar- 
row band imaging, a selection mechanism which is unbi- 
ased in colour-magnitude space. Spectroscopic observations 
were taken as follow up. The peri odic variable s (Tab le 9) 
came from lLamm et all \2004 ) and lLamm et all (|2005l ). We 
have used the entire perio dic variable c a talogu e, whereas the 
final membership list of iLamm et al.l (|2005l ) imposed fur- 
ther isochrone driven selection. Our photometry has been 
supplemented with the pho tometry for bright stars from 
iMendoza V. fc Gomed (|l980), as these stars were saturated 
in our catalogue. We include only those stars with probabil- 
ity of mem bership greater than 90 perce nt from the proper 
motions of (IMendoza V. fc Gomej|l980T ). We have also ex- 
cluded W33; the anomalous pho t ometr y of this star was also 
noted in IMendoza V. fc Gome3 i|l980l ). The full CMD with 
each form of membership criterion added is shown here as 
Figure IAT1 



A2 NGC2547 



A5 a Ori 



Ijeffries fc Oliveiral (|2005l ) carried out spectroscopy (Table 
5) of a photometrically selected sample. The photometric 
sequence was clear of the contamination in this instance and 
so the colour-magnitude selection should not have produced 
a bias in the CMD space. Radial velocity (RV) was used to 
define cluster memberships. The presence of Li below the Li 
boun dary was also used. Th e X-ray members (Table 4) are 
from I Jeffries fc Tollevl (|l998T ). 



Spect rosco pic members were taken from | Kenvon et al.l 
(2005) and lBurningham et al l (|2005bh . IKenvon et al.l (|2005l ) 



obtained spectroscopy of a photometrically selected sam- 
ple using a "close" selection around the suspected sequence 
in a CMD, which is potentially photometrically biased. 
Memberships were then confirmed on the basis of measure- 
ments of RV and EWs (Equi valent Widths) of Nal con - 
sistent with the group mean. iBurningham et all l|2005bl) 
however, used a "broad" colour-magnitude selection. The 
memberships were also confirmed here via RV and Nal 
and Lil EWs. IBurningham et all l|2005bl ) showed that us- 
ing a "broad" selection does not reveal sign ificantly more 
memb ers, implying that using members from iKenvon et all 
(2005) doe s not significantly bia s our results. Of the mem- 
bers from IKenvon et alj l|2005h . only t hose satisfying all 
the crit eria were included. In the case of IBurningham et alj 
(2005b]) a selection of > 80% probability of member- 
ship was applied to ensure a high percentage of genuine 
group members (Ta b le 21) . X-ray positions were taken from 
(2004J). We included those for which 
( 2004) (Table 20) found no optical coun- 



Sanz-Forcada et al 



Sanz-Forcada et al 



terpart and then cross-correlated the p ositions with our cat- 
alogue. The ISanz-Forcada et al.1 |200J) positions were used 
in preference to ROSAT d ata due to the highe r spatial res- 
olution of XMM-Newton. Ijeffries et ail <|2006f > shows that 
PMS members of this association are separated into two 
kinematic subgroups by heliocentric radial velocities, hav- 
ing different ages. In this work all but four stars from the 
sequence members are at a declination of less than 5=-2 
18' 00.0". This area is dominated by members from group 
2 as defined in Ijeffries et al.l i|2006l ). The resulting CMD is 
displayed as Figure [ 
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CepOB3b 



Figure A2. The full catalogue for the a Ori field (dots). Circles 
are members from lBu rningham et al ] d2005bt) , while asterisks are 
X-ray sou rces from | Sa rrz^Fbre ada^t**al.l J2004T) . Crosses are mem- 
bers from lKenvon et al .1 (|2005h . 



A6 Cep OB3b 

Spectroscopic members were taken from lPozzo et all (|2003T ) . 
They used a broad colour-magnitude selection with sub- 
sequ ent spectros c opic measurements of Lil EW, Ha and 
RV. H ozzo et alj (|2003h then used the spectroscopic mea- 
surements in conjunction with X-ray measurements to com- 
pare each star to the group mean and assign memberships 
(Tjible 14). Ha memberships were taken from lOgura et alj 
(2002) (Table 17), who used a photometrically unbiased 
narro w band imaging s urvey . The X-ray data were taken 
from INavlor fc Fabianl fl999), supplemented with the sec- 
ond ROSAT PSPC catal ogue. The latter conta ined two 
pointings; the one from INavlor fc Fabianl (|l999l ). and a 
second, non-over l appin g field. Both the unpublished and 
INavlor fc Fabianl l| 19991) catalogues were therefore used (Ta- 
ble 14), but where poss ible we took the positions from the 
INavlor fc Fabianl (I1999T) re d uctio n. In addition the X-ray 
sources from iGetman et al] ([2006) have been included (Ta- 
ble 18). These sources are from the CHANDRA ACIS detec- 
tor, and subsequently have a much higher spatial resolution 
and sensitivity than ROSAT. Both the ROSAT and CHAN- 
DRA data have been used as the latter is centralised on a 
smaller field of view about the cluster core, while many of 
the ROS AT detections are ou tside this area. In the cases 
where a IGetman" et al] (2006) source lies within the posi- 
tional error box of the sources from the ROSAT catalogues, 
the CHANDRA ACIS positions have been used. Periodic 
variables come from Littlefair et al (in preparation) (Table 
16). 

A7 IC348 



Spectroscopic and Ha data were taken from IHerbid |l998). 
He used photometry and a wide field Ha survey (grism spec- 
trograph), discarding Ha EWs < 2A(Table 27). Follow-up 
spectro s copy was undertaken on 80 stars within the field. 
IHerbid (1 19981 ) classified the spectra by comparison with 
dwarf spectral standards, using these to derive extinctions 
(Table 26) . The second set of s pectroscopic mem bers were 
taken from lLuhman et al] l|2003l ) (supplemented bv lLuhmanl 




Figure A3. The full 



catalogue 



f or the Cep OB3b fiel d 



(dots). Asterisks are X-ray sources from INavlor fc F abian (1999) 
IGetman et~aH J2006t) and th e Second ROSAT PSPC catalogue. 
Triangles are members from IPozzo et al ] l|2003l l. Circles are Ha 
sources from lOgura et al. I l l2002h . Crosses are the periodic vari- 
ables from Littlefair et al (in preparation). 



ll 19991). iLuhman et all (|2005h . iLuhman et aD (l2005h and 



ILuhman et al] (|2003l )) (Table 26). In ILuhman et al] ( 2003) 



spectral types were assigned using various models and then 
extinctions calculated by deredenning the st ars onto stan - 
dar d colours. The individu al extinctions from IHerbid (Il998l ) 
and ILuhman et aL (120031) are used he r e. X-r ay members 
were taken from IPreibisch fc Zinneckerl (|2002l ) (Table 24), 
where Ha and CHANDRA data were p r esente d with extinc- 
tions included. IPreibisch fc Zinneckerl (|2002l ) also include 
extinctions from a range of l iterature sources. We s upple - 
mented the X-ray data from IPreibisch fc Zinneckerl (|2002l ) 
with a portion of the Second ROSAT PS PC catalogue (Table 
25). Periodic variab l es wer e taken from I Cohen et aD (2004) 
and iLittlefair et al] (|2005l ) (Table 23), both being derived 
from / band wide-field surveys. The sequence with the mem- 
bers is displayed as Figure lA4l 



A8 A Ori 

Likely PMS members were selected from an R, R-I CMD. 
This sample was then spectroscopically observed using Li 
EW as an indicator of membership. The have used only those 
stars showing significant evidence for Li. 



A9 NGC2362 

1 J 1 

Dahml (|2005l ) used a wide field Ha survey to identify a sam- 
ple of 200 stars, which lay above the ZAMS in an optical 
CMD. This sample is then used for follow up spectroscopy, 
using Ha EW and Li EW as membership criteria. The mem- 
bers so identified are used here. 



A10 IC5146 

iHerbig fc Dahml ([20021 ) used a wide field survey to iden- 
tify any stars showing W(Ha) > 5 A. We used stars with 
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B2 NGC2547 
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Figure A4. The full catalogu e for the IC348 field (dots) Th e as- 
terisks are X-ray sources from IPreibisch fc Zinneckerl l|2002l l and 
the Second ROS AT PSPC cat alog ue. Circles are the peri odic vari- 
ables from lCohen et al.l lliooih andlLittlefair et al.l l|2005l ). Crosses 
are Ha sources from iHerbief ll 19981). Triangl es are s pectr oscopic 
memb ers with extinctions from lLuhman et al ] l |2003h and [Her big 
jl998h . 



IClarial f|l982h undertook a photoelectric survey in UBV and 
found 22 probable radial velocity members with (B-V) < 
0.1. Using individual reddenings to MS intrinsic colours he 
derived A v = 0.1 92, from E(B-V ) = 0.0 6 ±0.02 mags, with 
d m = 8.27 mags. iRobichon etail (I1999T ) used HIPPARCOS 
parallax data to calculate a d m = 8.18^q 26 P c - This later 
distance modulus used an alternative method to isochronal 
or intrinsic colour fitting, so will be most useful here. The 
studies used for memberships in Section IA2I have all con- 
firmed within uncertainties the parameters found in these 
two cases. The carried forward values are; d m — 8.18lQ 2g 
mags and Ay = 0.192. 



B3 ONC 



H2O MASER measurements are used in iGenzel et al1l|l98lf ). 
where expansion velocities of molecular clouds were used to 
infer the distance modulus. A result of d m = 8.3 8±0.37 mags 
was ob tained. This along with extinctions from iHillenbrandl 
<|l097f ) were used for this field. 



B4 NGC2264 



known spectral types and exhibiting a W(Ha) > 5A, dis- 
carding any sta r s lyin g below the Pleiades MS following 
iHerbig fc Dahinl <|2002h . 



All NGC6530 



X-ray members are from IPrisinzano et al] (|2005h using the 
CHANDRA ACI S detector. Ha members are taken from 
ISung et al l i|2000r i. 



APPENDIX B: LITERATURE DISTANCES 
AND EXTINCTIONS 

When obtaining parameters from the literature, as far as 
possible we avoided using values derived from PMS (Pre- 
Main Sequence) fitting, since we wish to derive ages inde- 
pendent of existing PMS isochrones. However, values derived 
from MS fitting (i.e fitting the upper part of the sequence) 
were used as these are independent of assumptions about 
the age. The following section details the literature sources 
used for each field. 



Bl h and \ Per 

Since the catalogues were combined and normalised to % P er i 
only the parameters for this cluster we re required. The pa- 
ramet ers for this cluster were taken from lCapilla fc Fabregati 
42002} ), A v = 1.71, d m = 11.7 ± 0.1 mags. This study fitted 
narrow-band photometry with MS isochrones. The values 
of extinction and distance mod ulus are in line with other 
recent determinations such as: iMarco fc Bernabeul (l200ll) 
E(b- y) = 0.39 ±0.05 and d m = 11.56 ± 0.20; ISlesnick et all 
(2002) d m = 11.85 ± 0.05 and E(B-V) = 0.56 ± 0.01. 



In iDahm fc Simonl (|2005l ) a summary of values from pre- 
vious studies is presented along with the methods of cal- 
culation. All the entries listed are from forms of isochrone 
fitting, bu t we consid e r only those studies fitting the main 
sequence. ISung et all (|l997h used colour-magnitude selec- 
tion in R — Ha, then fitted the ZAMS to assumed B type 
stars. An average distanc e modulus o f 9.4 ± 0.25 mags 
and A v = 0.23 resulted. iPerez et all (|l987h used shifts 
of individual OB type stars to intrinsic colours leading to 
A v = 0.19 and d m = 9.88 ± 0.17 mags. iPark et al.1 (|2000D 
selected stars using Ha and X-ray emission, then used MS 
isochrone fitting with the bright stars, deriving a distance 
modulus of 9.4 ± 0.10 mags and E(B-V) = 0.066 ± 0.034. 
An average of these values has been taken being d m — 9.6 
mags and Ay = 0.21. 



B5 a Ori 



iBrown et al.l l| 19941 ) found parameters using the lower part 
of the sequence, without the use of isochrones. Derivations 
of surface gravity from narrow band photometry were cali- 
brated to intrinsic colours of the theo r etical and empirica l 
models in IStraizvs fc Kurilienel l}l98ll ). IBrown et all (|l994l ) 
used MK spectral types in conjunction with effective tem- 
perature and absolute bolometric ma gnitudes to calculate 
surfac e gravities and radii of stars. In IStraizvs fc Kurilienel 
l|l98ll ) however the calibrations were done without taking 
into account PMS or post MS objects. IBrown et al.l (Il994l ) 
applied these calibrations to all spectral types, clearly at the 
age of around 3 Myr many of these stars are indeed PMS 
objects. A set of subgroups were selected and the subgroup 
lb in this paper covers our field. The values for this sub- 
group were given as d m = 7.8 ± 0.39 mags, with a value 
of E(B-V) = 0.06. The distance modulus from HIPPAR- 
COS parallax measurements is quoted as d m = 7.73l"g'gQ 
to the star a Orionis. The large uncertainties make this 
value unusable. Therefore the values taken forward were 
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from iBrown et all (|l994h . namely d m = 7.8 ± 0.39 mags 
and E(B-V) = 0.06 (Ay = 0.18). 

B6 NGC7160 

ISicilia-Aguilar et al.1 (120041 ) and ISicilia-Aguilar et ail (|2005l ) 
used membership criteria described in Section IA4I These 
papers used spectroscopically derived spectral types and in- 
trinsic colours to calculate the extinctions, and MS isochrone 
fitting for the distance modulus. The resulting values were 
A v = 1.17 and d m = 9.77 mags. 

B7 Cep OB3b 

iMoreno-Corral et alJ (| 19931 ) used optical and IR photometry 
of O, B and A type stars. These were compared with MS 
calibrations to calculate individual e xtinctions and dist ance 
moduli, which were then averaged. iPozzo et al.1 (120031) re- 
calculated these parameters from the Moreno-Corral et al.l 
1 19931 ) individual ex tinctions and distan ce moduli, but ex- 
cluded star 11 from iBlaauw et all (|l959l ) which is of dubi- 
ous membership. The resulting values which we also adopted 
were d m = 9.65 ± 0.20 and A v = 2.81 ± 0.10. 



used, with the average extinction being applied to the re- 
maining stars. Due to the need for individual extinctions, 
as with the ONC and IC348 the extinctions were applied 
before fitting. 

B12 NGC6530 

IPrisinzano et al] (|2005l ) used MS fitting of the stars in the 
blue envelope. These are MS stars assumed to be at the clus- 
ter distance or closer. They c alculate d m ~ 10.4 8 and use the 
value of the extinction from ISung et al.1 l|2000l ). ISung" et al.l 
(|2000h fitted a ZAMS relation individually to 30 early-type 
stars, yielding two different average distance moduli, assum- 
ing one was due to binarity they adopted the larger dis- 
tance modulus, d m ~ 11.25 ± 0.1 mags. The Q-method is 
then used for these ear ly- type stars yielding Ay ~ 1.09. In 
IPrisinzano et al] (|2005l ) a useful summary table of the dis- 
tance moduli calculated to this cluster is presented along 
with a summary of the literature of the field. For this study 
we have used both d m ~ 10.48, d m ^ 11.25 and used 
A v ~ 1.09. 



B8 IC348 

We adopted the distance modulus of IStrom et al.1 l|l974h . 
who used UBVKL photometry of 20 stars with V < 14 be- 
lieved to be cluster members, in conjunction with various 
literature spectral types and MK classifications to derive a 
distance modulus of 7.5lg'i6 ma g s - This is the same dis- 
tance modulus as that derived to the Per OB2 association, 
d m — 7.5 mag s. A ful l discu ssion of the distance to IC348 
is presented in iHerbid (|l998l ). The extinction s for th i s field 
were individual star extinctions and were from lHerbid (Il998l ) 
and iLuhman et al] i|2003h . 



B9 A Ori 

We have used the HIPPARCOS distance modulus to this 
cluster, dm = 7.90±0.17. The extin ction to this cluster is cal- 
culated in iDiplas fc Savage! i| 19941 ) as E(B-V)=0.12, equat- 
ing to Av — 0.36. 



BIO NGC2362 

iBalona fc Lanevl (Il996l ). using narrow band photometry of 
B-type stars d e rive A y = 0.31 and d m = 10.87 ± 0.03. 
iMoitinho et al.l (|200ll ) also used the B-type stars to obtain 
Ay — 0.31 and d m = 10.85. The values adopted here are 
10.87 ±0.03 and Ay = 0.31. 



Bll IC5146 

iHerbig fc Dahinl (2002) derived a distance modulus from 
a comparison of the early-type stars with models yielding 
dm — 10.4. They then calculated the reddening by assuming 
their distance modulus and using the known spectral types 
of a selection of 38 stars, to yield Ay — 3.0 ± 0.2. We have 
adopted d m — 10.4. The individual extinctions used for the 
B type stars to calculate the distance modulus have been 



Field & CCD (field.ccd) ID RA 



DEC 



x position (CCD) y position (CCD) MAG UNCERTAINTY FLG COL UNCERTAINTY FLG 



1.04 



173 02 22 19.173 +57 08 11.19 978.686 



2770.194 



12.230 0.010 



OO 



0.386 0.016 



OO 



Table Bl. The members for each field are available on the cluster home page and in the electronic version of this paper; this table is a sample. 



Field 


Ay 


Source 


Method 


d m (mags) 


Source 


Method 


Age(Myrs) 


Source 


NGC2547 


0.15 


1 


Q-method 


o -,0+0.29 
S - i8 -0.26 


11 


HIPPARCOS 


30-45 


27 


X P er 


1.57 ±0.079 


2 


Q-method 


11.7 ±0.1 


10 


MS fitting 


12.6 


2 


NGC2264 


0.371 


4 


Q-method 


9.6 


13 


Avg of many values 


3 




NGC7160 


1.17 


5 


Intrinsic colours 


9.77 


14 


Isochronc fitting 


10 


5 


Cep OB3b 


2.882 


6 


Q-method 


9.65 ± 0.20 


15 


Spectroscopy 


5.5 


20 


a Ori 


0.18 


8,9 


Iso fitting 


7.8 ±0.39 


17 


HIPPARCOS 


3 


17 


ONC 


A* 


19 


Spectral classes 


8.38 ± 0.37 


18 


H 2 MASER 


« 0.8 


19 


IC348 


A* 


7 


Spectral classes classification 


7 c+0.14 
' -°-0.16 


16 


HIPPARCOS 


2-3 


7 


A Ori 


0.37 


24 


Spectral classes 


7.90 ±0.17 


24 


HIPPARCOS 


4, 6-7 


28 


NGC2362 


0.31 


25 


Average of literature values 


10.87 ±0.03 


26 


MS fitting 


~ 1.8, 3.5-5 


25 


IC5146 


3.0 ±0.2 


23 


Spectral type classification 


10.4 


23 


MS fitting 


1 


23 


NGC6530 


1.09 


22 


MS fitting 


10.48 


22 


MS fitting 


2.3, 1.5 or 0.5-1.5 


29, 30 and 31 



Table B2. Extinction, d m , ages and sources. A* = Individual values from sources indicated. 



Field 


Q-method 


Literature 


X and h Per 


1.571 and 1.728 


1.82 and 1.71 


NGC2264 


0.371 


0.22 


NGC2547 


0.15 


0.192 


Cep OB3b 


2.882 


2.81 ±0.1 



Table B3. Extinctions (mags) from Q-method and from the literature. 



1 IClarial l|l982h. 2lslesnick et al.1 J2002h. 3lSanner et all d2000h. 4 fMendoza V. fc Gomej h980h 5 ISicilia- Aguilar et all (1200411. 6 Ijordi et alj lll996r i . 7lHerbid lll99Sft andlLuhman et al.l 

f ), 8lBeiar et all <199S|) . 9 lSherrv et all feOMhlBrown et alj Jl994h . lOlCapilla fe Fabregaj <2002i). lllRobichon et all Jl999h, 12lSanner et al.l feOOd). 13 Dahm fc Simonl ll2005h, 14 
i!-Aguilar et alj Il2004) 15lPozzo et alj koOci). 16[Herbiel lll99slvi7lSherrv et al.1 J2004h. 18lGenzel et alj lll98lJ). 19 Iffillenbrandl <199^. 20ljordi et all ll996h . 2l|PipIas fc Savage 
). 22 |Prisinzano et al.1 ll2005h, 23lHerbig fc Dahml hoO^ i, 24 lDiplas fc Savage! jl994l) , 25 |Dahml ll2005h, 26lBalona fc Lanevl jl996h , 27 fjeffries fc Oliveiral ll2005f r28 lMurdin fc Penston 
) and lDolan fc Mathieul (|200lh . 29 ISiess et alTl|2000l ), 30 Isung et al.l feOOO ) "and 31 Ipamiani et al.l [|2004h The used values of extinction for the lower age sequences is largely 
inconsequential as the isochrones are parallel to the dereddening vector at the lower mass end. 
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